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Fully relativistic first-principles electronic structure calculations based on a noncollinear local spin density 
approximation (LSDA) are performed for pyrochlore iridates ¥2^267 and Pr2lr207. The all-in, all-out antifer¬ 
romagnetic (AF) order is stablized by the on-site Coulomb repulsion U > Uc the LSDA-i-t/ scheme, with 
t/c ~ 11 eV and 1.3 eV for YfLciOi and Pr2lr207, respectively. AF semimetals with and without Weyl points 
and then a topologically trivial AF insulator successively appear with further increasing U. For U - 1.3 eV, 

¥2^267 is a topologically trivial narrow-gap AF insulator having an ordered local magnetic moment ~ 0.5//5/Ir, 
while Pr2lr207 is barely a paramagnetic semimetal with electron and hole concentrations of 0.016/Ir, in overall 
agreements with experiments. With decreasing oxygen position parameter v describing the trigonal compres¬ 
sion of Ir06 octahedra, Pr2lr207 is driven through a non-Fermi-liquid semimetal having only an isolated Fermi 
point of Fg, showing a quadratic band touching, to a Z2 topological insulator. 


Pyrochlore iridates A2lr207 [1] have attracted great interest 
for experimental observations of an itinerant chiral spin liq¬ 
uid [2] below 1.4 K and a possible underscreened Kondo effect 
around 20 K for A=Pr [3], a finite-temperature transition from 
a semimetal to an antiferromagnetic (AF) insulator for A=Nd, 
Sm, and Eu and a Mott insulator for a shorter ionic radius of 
A = Gd, • • •, Lu and Y [4-6]. Ir 5d electrons are also ad¬ 
vantageous for achieving topologically nontrivial states with 
and without electron correlations [7], because of a large rela¬ 
tivistic spin-orbit coupling comparable to an electron transfer 
and the onsite Coulomb interaction U, as experimentally ver¬ 
ified in Sr2lr04 [8]. Theoretical proposals include a Z2 topo¬ 
logical band insulator (TBI) [9-13] and Mott insulator [11], 
a Weyl semimetal and an axion insulator [14, 15] accompa¬ 
nied by all-in, all-out magnetic order [16, 17], a rhombohedral 
Z2 TBI under the [111] pressure [18], and a non-Fermi-liquid 
semimetal [19] associated with a quadratic band touching of 
electronic dispersions at the Fermi level [12, 20]. However, 
understandings of the underlying detailed electronic structures 
from the first principles remain open to date. 

One of the obstacles is that reliably treating localized rare- 
earth 4/ electrons at the A site is unfeasible in first-principles 
calculations. Calculations on ¥2^267 using crystal parame¬ 
ters for A = Pr, Nd, Sm, and Eu on the less correlated side 
may yield a discrepancy, because the ionic radius and thus the 
lattice constant a is longer by 1-3% in these materials than in 
¥2^267. It is also crucial to understand the dependence of the 
electronic structure on the only two crystal structure parame¬ 
ters allowed for A2lr207, i.e., the lattice constant a and the 
oxygen position parameter x [21], which vary with changing 
A site elements. It has been argued that even in the least cor¬ 
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related material Pr2lr207, the all-in, all-out AF metal (AFM) 
is stable within the local spin-density approximation (LSDA) 
based on the linear muffin-tin orbital method (LMTO) even for 
U = 0 [14]. This sharply contrasts to the experimental finding 
that Pr2lr207 is a paramagnetic semimetal (PSM) [2-5]. 

In this Letter, we provide extensive fully relativistic LSDA 
and LSDA-rf/ electronic structure calculations on ¥2^267 
(x=0.335, a =10.106 A) [22] and on La2lr207 with varying 
crystal parameters around the experimentally observed values 
(x=0.330, (2=10.400 A) [23] for Pr2lr207. A transition from a 
PSM to an all-in, all-out AF semimetal (AFSM) occurs at f/ = 
Uc with Uc ~ 1.1 eV and 1.3 eV for ¥2^267 and Pr2lr207, 
respectively (Fig. 1), in contrast to the previous LMTO results 
mentioned above [14]. 24 Weyl points [14] are observed at 
U = 1.15-1.2 eV and 1.35-1.5 eV for ¥2^267 and Pr2lr207, 
respectively. With further increasing U, the AFSM turns into 
a topologically trivial AF insulator (API). The API state of 
¥2^267 with U = 1.3 eV has an ordered local magnetic mo¬ 
ment mir ~ 0.5yU5/Ir, as shown in Fig. 1 (b), in reasonable 
agreement with the neutron-scattering experiments. [22] (The 
orbital and spin components are given by 0.29pB and 0.2 lyu^, 
respectively.) This value of U is smaller than U = 2.0 eV [24] 
in Sr2lr04, probably because the localized nature of Ir 5d elec¬ 
trons is pronounced by the layered structure in Sr2lr04. With 
the same value of U ~ 1.3 eV, Pr2lr207 is located on the verge 
of the PM phase. (See Fig. 1 (a).) The PSM state for Pr2lr207 
has a cuboidal electron Fermi surface centered at the F point 
and toric hole Fermi surfaces around the L points. With de¬ 
creasing X and thus increasing Ir-O-Ir bond angle from 121° 
(x = 0.35) to 148° (x = 0.30), we observe a transition to 
a Fermi-point semimetal hosting a quadratic band touching 
of two doubly degenerate bands which form Fg at the Fermi 
level. A further decrease in x drives a level crossing of Fg and 
Fg, resulting in a Z2 topological insulator for x ~ 0.31. 

The calculations have been performed by means of a fully- 
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FIG. 1: (color online) Phase diagrams of (a) Y 2 lr 207 and (b) Pr 2 lr 207 as functions of U, obtained with the experimentally observed crystal 
parameters; the total energy difference AS between paramagnetic semimetal (PSM) and AF semimetal (AFSM) or insulator (AFI) in the left 
axis, and the magnitude mir of the ordered local magnetic moment per Ir site in the right axes. 




FIG. 2: (color online) Electronic band dispersions (left) and atom projected densities of states (right) of (a) the all-in, all-out AFI for ¥ 2^267 
and of (b) the PSM for La 2 lr 207 with the crystal parameters of Pr 2 lr 207 , obtained by the LSDA-ht/ (U = 1.3 eV) and the LSDA (U = 0), 
respectively. Red, green, magenta and blue curves represent Y or La, Ir, 01 (48/), and 02 (SZ?) states, respectively. The Fermi level is denoted 
by grey dashed lines. 


relativistic non-collinear spin density functional theory within 
the LSDA and LSDA-l- U methods implemented in OpenMX 
code [25], which is based on the norm-conserving pseudopo¬ 
tential and a linear combination of multiple pseudo-atomic 
orbitals. (For computational details, see Supplemental Ma¬ 
terials.) We use (9,9,9) to (18,18,18) uniform ^-point mesh. 
The pseudopotentials and pseudo-atomic orbitals are severely 
checked by comparison with the all-electron full-potential lin¬ 
earized augmented plane-wave (FLAPW) method [26], so that 
our fully relativistic scheme is well compared to the FLAPW 
results in the LDA with the spin-orbit interaction being treated 
in the second variation scheme. (See Fig. SI in Supplemental 
Materials.) The total-energy calculation using OpenMX re¬ 
veals a structural stability around a = 10.4 A and x = 0.33 
for La 2 lr 207 (see Supplemental Materials), which are close to 
the experimental observations for Pr 2 lr 207 and justifies the La 
replacement as long as the low-temperature rare-earth mag¬ 
netism is irrelevant above 1.5 K[2]. Thus, we henceforth re¬ 
place Pr with La for calculations on Pr 2 lr 207 . We assume 
that the translational invariance is fully preserved, while the 
ordered magnetic moment directions have been fully relaxed 
during the calcuations without adding any energy barrier and 


converged to the all-in, all-out AF structure in all the magnet¬ 
ically ordered cases shown below. 

Figure 2 shows the band structure and atom projected den¬ 
sity of states (DOS) of (a) a narrow-gap AFI obtained for 
Y 2 lr 207 with the LSDA-hf/ (U = 1.3 eV) and (b) a PSM 
obtained for Pr 2 lr 207 with the LSDA (U = 0). The Fermi 
level has been chosen to be E’f = 0. In both cases, the bands 
in E e [-2,0.5] eV shown in the atom projected DOS (Fig.2) 
are mainly composed of anti-bonding states of Ir 5d and 01 
(48/) 2p electrons in Ir06 octahedra. The energy windows 
E € [0.5,2.0] eV for Y 2 lr 207 and E e [0.5,1.5] eV for 
Pr 2 lr 207 are given by a crystal-field gap due to the Ir06 octa¬ 
hedral coordination. 

In this AFI for Y 2 lr 207 , the charge gap is suppressed by a 
proximity to the AFWSM. Figure 3 shows an evolution of the 
electron dispersions around a Weyl point with increasing U 
from 1.2 to 1.3 eV. Weyl points are confined to the Y-X-L and 
symmetry-related planes [14] by the magnetic space group 
symmetry Ed3m' if they exist. They indeed appear at k = 
^(0.44,0.44,0.35) and (0.46,0.46,0.41) for C/= 1.15 eV and 
1.2 eV, respectively, and symmetry-related points, with their 
energy level being located about 35 and 24 meV higher than 
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the Fermi level, leaving electron and hole Fermi surfaces as¬ 
sociated with the concentration ~ 0.008/Ir = 1.3 x 10^^ /cm^. 
The Weyl point shifts towards the L point with increasing 
U and eventually disappear at t/ ~ 1.25 eV. This leaves 
an topologically trivial AFSM with similar electron and hole 
Fermi pockets. Then, an indirect gap ~ 5 meV opens be¬ 
tween the valence band top at the L point and the conduc¬ 
tion band bottom at the F point for U = 1.3 eV. It is likely 
that disorder effects localize thermally excited dilute carriers 
in this narrow-gap AFI and inelastic scattering events beyond 
the LSDA-rf/ shcheme give rise to the variable range hop¬ 
ping, as experimentally observed in Y 2 lr 207 [27], as well as 
in Nd 2 lr 207 [28] and Eu2lr207[29]. 

The PSM for Pr 2 lr 207 contains a twofold degenerate 
cuboidal electron Fermi surface around the F point and 
twofold degenerate hole pockets around the L points, as 
shown in Fig. 4. From the Fermi volumes, both the elec¬ 
tron and hole carrier numbers are given hy He = Uh = n = 
0.016/Ir = 2.22 x 10^^ jcvc? for Pr 2 lr 207 . From the Hall co¬ 
efficient Rh = -2.4 X 10“^ cm^C measured in Pr 2 lr 207 [3], 




the ratio of the electron and hole mobilities jde and jik is esti¬ 
mated to be a rather large value, = 13, through 


Rh 


(neiil-nhul) 

eiriefie+nhfih)'^ 


with the electron charge e. 


(a) 



E{T+){eV) 
- 0 

- - -0.05 

- - - 0.1 




a (A) 


FIG. 3: (Color online) Electron band dispersions in the T-X-L plane 
for ¥ 2 ^ 267 , (a) in the AFSM with U = 1.2 eV, (b) around the 
AFSM-AFI boundary at t/ = 1.25 eV, and (c) in the AFI with 
U = 1.3 eV. 



FIG. 4: (color online) Fermi surfaces of La 2 lr 207 for the crystal 
structure of Pr 2 lr 207 . 


FIG. 5: (color online) Phase diagrams of La 2 lr 207 : (a) the energy 
level F(rp - Fp and (b) the level difference AE = E{r^) - FfFg), 
which changes the sign at the white solid line. ■: a 242 PSM, •; 
a 242 paramagnetic Fermi-point semimetal (PFPSM), ♦: a 422 Z 2 
TBI, and ▲: a 422 PSM. The numbers 242 and 422 in the above 
represent those of the degeneracies of the three levels at the F point 
around Fp in the ascending order in energy. 

Now we study the dependence of the electronic structures 
on a and x within a PSM or a band insulator, as (a,x) in 
Pr 2 lr 207 depends on the temperature, a hydrostatic pressure, 
and possibly the sample quality. Figure 5 shows the phase 
diagram of La 2 lr 207 obtained with f/ = 0. It contains four 
phases; (a) a PSM as obtained for Pr 2 lr 207 , (b) a paramagnetic 
Fermi-point semimetal (PFPSM) of two doubly degenerate 
quadratic bands touching at the F point at F’p from above and 
below, (c) a Z 2 TBI, and (d) another PSM with even smaller 
electron and hole Fermi pockets than in Pr 2 lr 207 . This evo¬ 
lution of the electronic structure is attributed to a change in 
the single-orbital/Ir Kramers degenerate states closest to F’f, 
which form a closed manifold consisting of four doubly de¬ 
generate bands, as shown in Figs. 6 (a) to (d), respectively. In 
particular, the key role is played by an energy level crossing of 
a twofold degenerate F^ and a fourfold degenerate Fg, as we 
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show below. (The degeneracy of the latter is protected by the 
cubic and time-reversal symmetry, while that of the former is 
only by the time-reversal.) 

When X > Vc ~ 0.31, namely, in the upper part than a white 
solid line in the phase diagram (Fig. 5 (b)), the energy levels 
at the r point satisfy the relation £’(rp < E(r'^) < E(r^), as 
shown in Fig. 6(a). With decreasing x, £’(Fg) becomes shal¬ 
low below £’f, as seen from the color map of Fig. 5 (a). Ac¬ 
cordingly, both electron and hole Fermi volumes shown for 
a = 10.40 A and v = 0.33 in Fig. 4 gradually decrease. This 
slight increase of £’(Fg) appears to be due to a rapid increase 
of £’(Fp approaching £’(Fg) and £’f, as shown in the color 
map of Fig. 5 (b). Then, around v ~ 0.32, the electron Fermi 
surface eventually shrinks into a single point, while the hole 
Fermi surfaces disappear, yielding a PSM hosting a single 
Fermi point with £’(Fg) = £’f, as shown in Fig. 6 (b). This 
paramagnetic Fermi-point semimetal (PFPSM) appears in a 
white region in Fig. 5 (a). (To be precise, we obtained this 
solution for 9.8 A < <2 < 10.5 A in the case of v = 0.32 and 
for 10.0 A< a < 10.5 A in the case of v = 0.31, as marked by 
green points in Fig. 5.) 

On the other hand, with decreasing x below £^(rp be¬ 
comes larger than £’(Fg) in the lower part of the phase diagram 
than the white solid line (Fig. 5 (b)). This yields a narrow-gap 
insulator around v ~ 0.31. (To be precise, the insulating so¬ 
lution is found for x = 0.31 in the case of a = 9.8 and 9.9 A 
and for x = 0.309 in the case of a = 10.4 A, as shown in 
Fig. 6(c).) This is actually a strong Z 2 TBI with the Z 2 topo¬ 
logical indices [9] (1; 0,0,0), which has been confirmed by 
directly calculating the parity eigenvalues of all the occupied 
states at the time-reversal invariant momenta. The same topo¬ 
logical indices can also be obtained within the single-orbital 
manifold per Ir near £'f. Within this manifold, the occupied 
states at the eight time-reversal invariant momenta take the 
following irreducible representations; the fourfold degenerate 
Fg, fourfold degenerate -r at Z : (^, 0,0) and the two 
symmetry-related points, doubly degenerate (L| + Tp and 
at a single L point, e.g., (f, f, f), and (L~ L~) and L~ at 

the three symmetry-related points. Note that the Z 2 topologi¬ 
cal indices do not altered by any level crossing, for instance, 
of £’(Fp and £’(Fg), observed in the above phase diagram. 
Actually, the Z 2 TBI is realized by simply gapping out the 
band dispersions around the Fermi level, which certainly re¬ 
quires £’(Fg) < £’(rp, E(T^). With further decrease in x down 
to 0.30, we observe that bands cross the Fermi level, leaving 
small electron and hole Fermi pockets near the F and L points, 
respectively, for 9.8 A < a < 10.5 A, as shown in Fig 6(d). 
Note that the Z 2 TBI phase and the PFPSM phase may expand 
more if we employ a method which may improve the bandgap 
in LDA calculations, such as modified Becke-Johnson ex¬ 
change potential[30], hybrid functionals(e.g., HSE[31]), and 
GW methods[32]. 

The PSM solution obtained for the crystal parameters of 
Pr 2 lr 207 is clearly compatible with the experimental finding 
on the ordinary Hall constant [3], as we already discussed, 
if we assume that the mobility of the holes is an order of 
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FIG. 6: (color online) The (a, x) dependence of the band structure of 
La 2 lr 207 : (a) (a = 10.4 A, x = 0.33) the same structure as Pr 2 lr 207 , 
(b) (a = 10.4 A, X = 0.31), (c) (a = 9.8 A, x = 0.31), and (d) 
(a = 10.4 A, X = 0.30). 


magnitude smaller than that of the electrons. In this sce¬ 
nario, Shubnikov-de Haas oscillations with a single Lifshitz- 
Konsevich oscillatory component, experimentally observed in 
Pr 2 lr 207 [33], should be ascribed to the electron-like carriers. 
It will be possible to reach the quadratic band touching PF¬ 
PSM, if a compound with smaller x is found. This state can 
be experimentally verified by a V^-dependent Hall resistiv- 
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ity [19] and an angle-resolved photoemission spectroscopy. 
The bulk energy gap in our cubic Z 2 topological insulator 
phase can be as large as 20 meV, and is experimentally ob¬ 
servable, although it is required to find a material with even 
smaller v. 

When we almost finished the work, we noticed recent 
preprints [34, 35] on the LSDA-hDMFT calculations for 
A 2 lr 207 , which show a direct first-order phase transition be¬ 
tween the PM and an API in the space of U and the temper¬ 
ature. The stability of the paramagnetic state up to a mod¬ 
erate value of U is consistent with our results. However, Uc 
for the AF order is much larger than in our current results, 
which washes out the AFSM with or without Weyl points in 
the weak-to-intermediate range of U. The ordered local mag¬ 
netic moment, O.bSyU^ [34] or 0.51 jib [35], for Y 2 lr 207 and 
the charge gap, 0.3 eV for Y 2 lr 207 [34] or 0.4-0.5 eV for 
Nd 2 lr 207 [35], are also overestimated in comparison with our 
current results which agree with experiments. Note also that a 
recent angle-resolved photoemission spectroscopy measure¬ 
ment in the low-temperature all-in, all-out ordered state of 
Nd 2 lr 207 [16] has revealed half the symmetrized direct en¬ 
ergy gap of 50 meV at the L point [36], which is compati¬ 
ble with our results (Fig.3(c)) but not the LSDA-hDMFT re¬ 
sults [34, 35]. 
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Supplemental Materials 

COMPUTATIONAL METHODS 
OpenMX 

By using the OpenMX code[l], we have performed fully-relativistic first-principles band structure calculations based on the 
density functional theory (DFT) within the local density approximation (LDA)[2]. Fully relativistic pseudopotentials including 
spin-orbit coupling are generated by the Morrison-Bylander-Kleinman scheme[3], which includes a partial core correction[4]. 
The wave functions were expanded in a linear combination of multiple pseudo-atomic orbitals (LCPAO) generated by a con¬ 
finement scheme[5, 6 ]. Pseudoatomic orbitals are chosen to be LaS.O-^'SpSJl/l, Ir 7 . 0 - 5 ' 3 / 73 J 2 , and 05.Q-s?>p?>d\. The numbers 
after the chemical symbol (8.0, 7.0, 5.0) are the cutoff radii (bohr), while the latter parts of the expansion (e.g., s?>p?>d\f\ etc.) 
are the number of orbitals for the 5 ', p, d and / characters. We adopted the cutoff energy 250 Ry which was used in the calcula¬ 
tion of matrix elements associated with the difference charge Coulomb potential and the exchange-correlation potential, and in 
solving Poisson’s equation using fast Fourier transform (FFT). 


FLAPW 

We include spin-orbit coupling as the second variation after convergence of the scalar-relativistic self-consistent-field (SCF), 
i.e., we diagonalize the Hamiltonian matrix including spin-orbit interaction constructed from scalar-relativistic eigenstates and 
wavefunctions[10, 11]. We use Perdew-Zunger expression for exchange and correlation energy[12]. The k-space integration is 
achieved by the improved tetrahedron method[13]. For the face centered cubic cell, a (4,4,4) mesh is used during the SCF loops 
with 8 k point in the irreducible Brillouin zone (IBZ). Muffin-tin spheres radii assumed are 1.0 A for La, 1.1 A for Ir, and 0.8 A 
for O. Core states of 15* to Ad^^ for La, 5p^ for Ir, l^-^ for O fully relaxed by using the spherical part of the one-electron potential 
during the SCF iterations. The energy cutoff of plane waves used is 25 Ry. The charge density and potential are expanded using 
the cutoff of 100 Ry. 


CHECKING THE PSEUDOPOTENTIALS AND THE PSEUDO-ATOMIC ORBITALS 

In order to check the reliability of the pseudopotential, calculations were also done with the all-electron full-potential lin¬ 
earized augmented plane-wave (FLAPW) methods [7-9] with local density approximation (LDA) including spin-orbit inter¬ 
action. Figure SI shows comparison of the band structures of FLAPW methods (HiLAPW code) and OpenMX code, for (a) 
La 2 lr 207 with the crystal parameters of Pr 2 lr 207 and (b) ¥ 2^207 with the relativistic spin-orbit interaction. The calculated band 
structures obtained with the OpenMX code show good agreements with those with the HiLAPW code. This gives a support for 
the reliability of the present pseudopotential calculations and convergence in pseudo-atomic orbitals. 


STRUCTURAL STABILITY 

Figure S2 shows total-energy profile for La 2 lr 207 in terms of lattice constant a and internal parameter v for oxygen 48/ site. 
The calculated total-energy is lowest at the experimentally observed a and v for Pr 2 lr 207 . 
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FIG. SI: Calculated band structure of (a) La 2 lr 207 and (b) Y 2 lr 207 . Black lines are calculated by HiLAPW code. Red lines are calculated by 
OpenMX code. 
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FIG. S2: Total energy for La 2 lr 207 , the lattice constant a and internal parameter x for oxygen 48/ site changes 
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